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Abstract 

We examine the Type la supernova data in order to determine if it 
shows any signal of large scale anisotropy. The anisotropy is modelled 
by an extended Godel metric, which incorporates expansion along with 
rotation. The model is smoothly connected to the usual FRW type, 
while expressing anisotropic metric effects depending on certain pa- 
rameters. We find no significant signal of anisotropy in the data. We 
obtain bounds on an anisotropic redshift versus magnitude relation- 
ship, and accompanying parameters of the Godel-Obukhov metric. 



1 Introduction 



The large redshift type la Supernova data [Il[2l[3llll[5l[6l[7l[8] provides direct 
evidence for dark energy and has played a fundamental role in determination 
of cosmological parameters. In this paper we examine whether the data is 
consistent with the assumption that the universe is isotropic on large distance 
scales. There exist many observations which indicate a possible large scale 
anisotropy. These include the large scale anisotropy in the radio polarizations 
from radio galaxies [9l [10] , large scale alignment of quasar optical polariza- 
tions |TT], alignment of CMBR quadrupole and octopole [HI [13]. Here we use 
an anisotropic Godel type model to test for large scale anistropy. There is an 
interesting interplay between the redshift dependence and the anisotropy of 
the model which can be tested rather cleanly. We also test for anisotropy by 
using some simple extensions of the formula for distance modulus obtained 
in the standard Big Bang model. These extensions model the effect of a wide 
range of anisotropic metrics in the limit of weak anisotropy. In no case do 
we find a significant signal of anisotropy. Hence we use the supernova data 
to impose limits on these anisotropic models. 

The Godel universe [H] has received significant attention in the past few 
years, following the discovery that it was an exact solution of a string theory 
[I5l[THl[II|. Indeed since its discovery in 1949, the Godel universe attracted 
interest due to its fascinating properties, including explicit rotation, the ex- 
istence of closed time-like curves (CTCs), and an anti-Machian nature. It 
is well-known that Godel's original universe is an inappropriate model be- 
cause of absence of expansion and redshift. A number of Godel type metrics 
have been proposed which allow for expansion of the universe, and in which 
the quandary of CTC's are banished. In Ref. [TSl jTH] Obukhov generalizes 
Godel's model to include expansion, listing conditions for CTC's to exist. 
The revised model contains parameters which smoothly interpolate between 
the usual Friedman-Robertson- Walker (FRW)model and Godel's. The inde- 
pendent nature of vorticity and rotation, as discussed by Obukhov [HI [19] 
does not allow limits on rotation to be placed by limits on vorticity. Among 
distinctive predictions of the model are effects on the propagation of light. 

Comparisons with observational data have been scant. In 1982 Birch [9] 
claimed the existence of a large scale anisotropy in radio wave polarizations 
from extragalactic sources. He suggested that this might be observational 
evidence for universal rotation and anisotropy. Birch's controversial work 
was revisited by many authors with conflicting outcomes [20] • Studies dis- 
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missing the effect invariably use a statistic of even parity, while the effect 
itself is a corkscrew twisting or polarization with odd parity [TH]. In Ref. 
|10j . Jain and Ralston explored tests not requiring redshift information and 
found significant signals of anisotropy in a large sample of data. Several 
other observables of radiation propagating on cosmological scales have been 
found to indicate a preferred direction, all of which are aligned along the 
same axis [131 [121 HI]- In these studies no attempt was made to identify 
the physical origin of anisotropy. They may represent effects independent of 
gravitation and restricted to modifications of the electromagnetic sector, in 
which polarization observables are exquisitely sensitive. 

"Dark energy" is the popular motivation to consider models beyond the 
FRW metric. Current CMB data fits a fiat cosmology with f2jv/ ~ 0.27, 

~ 0.73 for the dark matter and dark energy density, respectively. These 
fits assume an isotropic universe. Meanwhile the alignment of low angu- 
lar momentum CMB multipoles, which naturally have directional features, 
actually contradict isotropy at a level greater than 2 sigma |12l [13]. This 
suggests we examine the data for Type la supernovas for consistency with 
isotropy. There is no model-independent way to go from effects on the CMB, 
and supernova effects in the current era, so that the two studies are inde- 
pendent. The specific Godel-Obukhov metric that we study does not predict 
any anisotropic effect on the CMB temperature due to the absence of shear 
in this model [TS] . 

It would be hard to take current supernova magnitude data and establish 
an effect of anisotropy without using a cosmological model. This is because 
the data comes from more than one group, and is necessarily affected by 
uneven sampling on the sky and in redshift. It is slightly different to ask 
whether the data supports isotropy. An invariant correlation (defined below) 
shows that supernova redshifts are strongly correlated with their angular 
positions. This is reasonably explained by selection. It turns out that su- 
pernova magnitudes are even more correlated with angular position than the 
correlation in redshift would suggest. Given the difficulty of establishing an 
anisotropic effect, we don't even attempt to assess the importance of this 
signal. The problem is that the relative importance of angular correlations 
of magnitudes is coupled to the model of redshift dependence, and vice- versa. 
We have not seen this discussed before: it might be that by incorporating 
anisotropic effects, the values of VLm or would change. Or, one might be 
led to place bogus hmits on anisotropy hinging upon particular assumptions 
about the redshift dependence. The resolution cannot be made on the basis 
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of a single test statistic and hinges on more detailed study we will describe. 
The outcome also depends on what are used for host galaxy extinctions in 
the supernova data, as we discuss. 

In the end our most reasonable fits do not show any signal requiring 
anisotropy. Yet somewhat surprisingly, the limits on anisotropy violation 
parameters are comparatively weak. That is, some small anisotropy cannot 
be ruled out. 



The Godel-Obukhov line element fl8[ [19] is 

ds"^ = dt^ - 2^f^)R{t)e'^'''dtdx2 - R\t){dxl + /te^^^Mx^ + dxl). (1) 
Let us comment on various terms: 

• Spatial axis X2 is singled out as breaking isotropy. In work below, we 



• The expansion factor R{t) is similar to what appears in the FRW met- 
ric. Since R(t) is a single scale factor for all spatial coordinates, it also 
appears in the coefficient of the off-diagonal term dtdx2- 

• Other parameters in the metric are a, k and m. Causal structure of 
space-time is regulated by the parameter k. For A; > there are no 
CTCs, and conversely k < has CTCs. 

• The classic Godel metric is obtained by setting R(t) = l,a(t) = 1, m = 
1, k = -1/2. 

• Parameter a determines the magnitude of acceleration of the fluid ele- 
ment due to rotation of the universe determined by the metric \W\ . 
Similarly, parameter m can be identified as proportional to vortic- 
ity. The relation of matter to the metric of course depends on the 
model, which has been examined in the context of Einstein's general 
relativity [21] |22| [23] as well as the Einstein-Cartan theory [19]. It is 
convenient to lump together the combination important for our study: 
the anisotropy parameter p defined by 



2 The Godel-Obukhov Model 



designate the axis of anisotropy more generally as A. 




(2) 
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2.1 Magnitude Versus Redshift 



The metric Eq. [T] predicts an anisotropic relationship for magnitude as a 
function of redshift for standard candles. Using the Kristian-Sachs [24J ex- 
pansion, Obukhov [TS] derives the following expression for the bolometric 
magnitude rriboi 

rriboi = rhhoi - 5 log^o (1 + p cos a) . (3) 

Here fhboi is the usual FRW bolometric magnitude and a is the angle between 
the source and the X2 axis. If we expand the second term on the right hand 
side to leading order in p it gives the "dipole" correction due to rotation that 
we will bound. In arriving at Eq. [3] we have kept only the leading order term 
in the redshift z which depends on the angular positions. This term is in fact 
independent of z. The term linear in z is found to be proportional to ujq/Hq, 
which is negligibly small; current limits imply < 10~^ [T8l[9]. Similarly, 
we ignore tiny kinematic effects from our reference frame's velocity relative 
to the rest frame of the CMB. If p << 1, which indeed must be the case, we 
may expand the logarithm and keep only the leading order term in p. 



2.2 Relation to Distance Moduli 

The distance modulus Pp of a supernova is defined hj Pp = m — M, where 
m and M are the apparent and absolute magnitudes respectively. This is 
related to the luminosity distance by the relationship 



Pp = 5 logio + 25 



(4) 



We parametrize the luminosity distance in the FRW limit, setting a = 0, as 
follows [25] 



l + ^(l-go)^-^(l 



(5) 



where go and jo are respectively the deceleration and the jerk parameters, 
defined in Ref . [23] . We can now express the distance modulus in the Godel- 
Obukhov metric, Eq. [H as follows: 



yUp = Ppo + 5 logio ^ + 5 logio 
- 5 logio [1 + P cos a] , 



1 1 



(6) 
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Figure 1: Scatter plot of angular positions of the 186 object gold and silver 
set compiled in Ref. [26] in Equatorial J2000 coordinate system. Symbols 
indicate low redshifts {z < 0.4, plusses) and high redshifts {z > 0.4, squares). 

where 

cos a = A • f . (7) 

Here f is the direction of the source, and (to repeat) A is the preferred axis 
of anisotropy of the metric. We have included an overall parameter /ipo to 
allow for a zero point correction, which includes uncertainty in the Hubble 
parameter. 

3 Angular Correlations 

Here we test for signals of angular dependence in the type la supernova 
data. We examine the "gold" data set comprising 157 sources from Ref. 
[26] . as well as the gold and "silver" data set of 186 sources. A scatter 
plot of the data is shown in Fig. [H One sees that the data population is 
quite unevenly distributed in the sky, particularly at high redshifts. There is 
nothing remarkable about this distribution, of course. 
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In many statistical tests the sample selection might cause a spurious signal 
of angular dependence. Hence naive signals of anisotropy must be interpreted 
with care. Conversely, selection effects may also mask a true signal of angular 
dependence that might be present. We now turn to tests of correlation that 
do not depend on the sample population. 

3.1 JM Correlation Statistics 

We quantify correlations of distance moduli with angular positions on the 
sky by using the Jupp-Mardia [27] (JM) angular correlation coefficient r. 
In general, one maps angular variables into "vectors" (upper case Xj) that 
transform linearly under a change of coordinates. From the map statistics 
which are rotationally invariant are constructed. 
The procedure uses correlation matrices 

CXX,^, =< {X - XUX - X), >■ CxY,^3 =< {X - X^Y - Y), >, 

etc. where < > indicates a summation over the data set and X, Y refer to 
the mean values over the sample. Explicitly < XiXj >= J2a X^Xj', where 
X" is the jth vector component of the a-th data point. The general JM 
statistic r^(X, Y) is the combination 

r\X, Y) = tri C^\CxYCy].CYx ), (8) 

where matrix multiplication is implied, denotes an inverse matrix, and 
tr indicates the trace. 

The JM statistics are a very general test of correlation of the underlying 
distribution /(X, Y). Distributions are uncorrelatedif 

/(X, F) = /(X)/(r), {null model), 

where /(X) is the marginal distribution of X, and similarly for Y. For a 
sample of size n, JM show that the product nr^, where is the statistic 
defined in Eq. |H1 from an uncorrelated null distribution is distributed like 
xl, where u is the total number of components of X and Y. Extensive Monte 
Carlo calculations confirmed this distribution in Ref. [TO] . 

It is very important that JM test statistic is independent of the partic- 
ular marginal distributions. Population (marginal distribution) effects are 
automatically removed by inverse matrices used in the definition. For our 
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gold and silver (z) 
gold and silver (/x) 
gold and silver {Ay) 



sample 
gold (z) 
gold (/x) 
gold (Av) 



157 
157 
133 
186 
186 
156 



n 



2 

nr 
20.4 
27.1 
15.5 
16.3 
23.5 
9.5 



P — value 
1.4 X 10-^ 
5.6 X 10-6 
1.4 X 10-3 
9.8 X 10"^ 
3.2 X 10-^ 
2.4 X 10-2 



Table 1: Jupp-Mardia invariant statistics for the correlation of redshifts 
(z), distance moduli (/i), and host extinctions (Ay) of supernovas with their 
angular position on the sky. Symbol n is the sample size, and the P-value 
is the probability to see the statistic or larger ones in an uncorrelated null 
distribution. Supernovas for which host extinction values are not available 
are removed while computing the statistics for Ay- 

purposes we assign f ^ X, where f is the unit vector of the position on the 
sky in equatorial coordinates. We assign scalar quantities such as distance 
moduli Hp Y. Eq. [8] applies without change, the meaning of C{YY) 
simply becoming the sum of squares of Y, and so on. By construction (and 
previous consistency tests) the uneven sampling of supernovas on the sky 
cannot create a significant false correlation of nr"^ if the sample distribution 
obeys the null model. Thus the JM test statistics do not require any partic- 
ular model of anisotropy, but they can rule out a null hypothesis of isotropy. 
The following facts are found: 

• Supernova Magnitudes Break Isotropy: Evaluation of the JM correla- 
tions produced Table [1] comparing nr'^{z) and nr'^i^jj). As discussed in 
the introduction we cannot expect supernova sampling to be uniform 
in "depth", namely redshift. This is reflected in the large correlation of 
redshift with angular position, which indicates sampling effects. Mean- 
while in both the gold set, and the gold and silver set, the magnitudes 
themselves are much more correlated than the redshifts. This cannot 
be attributed to sample skewing by the Milky Way: the numbers rep- 
resent a correlation of magnitudes with position, not the distribution 
of the positions themselves. Nor can the results be explained away 
by extinction effects of our galaxy, unless those effects have been es- 
timated incorrectly. Our own galaxy's extinctions have already been 
removed by the observers. We are not sure how to assess the relative 
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probability of the magnitude correlation given the redshift correlation. 
For reference, the ratio of P — values is of order 0.03 — 0.04, perhaps 
a significant ("2-cr") effect. 

• Host Galaxy Extinctions Break Isotropy: In Ref. evidence was 
found for evolution or bias in the host galaxy extinctions of supernova 
in the gold set. The trend is that host extinctions are strongly corre- 
lated with their deviation from the Hubble plot. The method by which 
extinctions are assigned by observers tends to produce the same signal 
as acceleration. We decided to compute the JM correlation of host ex- 
tinctions with angular positions (Tabled]). For the gold set the result 
i£l 

nr^{Av) = 15.5, 

a P-value of 1.4 x 10^^. This "3cT"-effect is statistically quite significant. 
Perhaps the combination of effects explains the pattern seen in Table 
[H Admittedly there is no reason a priori for the two effects to be 
oriented in a way to add. Nor can this possibility be checked with the 
JM correlations themselves. 

• Anisotropy of Magnitudes and Host Galaxy Extinctions Are Related: 
Comparing the correlation of host extinctions and magnitudes requires 
a directional statistic. We computed y-weighted average unit vectors 
< YX >, 

< >= ( X] 0^" " y)^" )normalized- 

a 

Again the sum is over galaxy labels a. Subscript -normalized indicates 
that result is normalized after making the sum. This statistic has two 
degrees of freedom, and like the JM correlations its direction is not 
strongly affected by the marginal distribution on the dome of the sky. 
By construction < YX > vanishes if distributions are uncorrelated, and 
any bias in the X marginal distribution should average symmetrically 
to zero in an uncorrelated null distribution. Moreover, the relationship 
of < YX > calculated for Y = z and for Y = Ay should be random 
if the host extinctions are not causing the extra angular correlation 

24 sources for which extinctions are not reported have been deleted. 
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of magnitudes seen in Table [TJ The naive probability of overlapping 
one unit vector with another isotropically is P(< YiX >, < Y2X > 
) = (1- < YiX > ■ < Y2X >)/2. We only conduct this study for the 
gold set and those entries for which Ay are listed. Surprisingly all the 
F-weighted vectors line up, with: 

P(< zX >, < fxX >) = 2.5 X 10-2 

P(< AyX >, <zX >) = 6.9 X 10-2 

P(< AyX >, < fxX >) = 1.2 X 10-2 

One of these statistics is dependent; as already discussed the z — fi 
coincidence is also perfectly expected. The direction of < AyX > is 
RA = 114°, DEC = 8°. It is hard to understand how the host extinc- 
tions might have become correlatecj^ in this way. However, given the 
fact that this angular correlation exists, the higher correlation of mag- 
nitudes reported relative to the correlation in z seems to be explained. 

Up-Down Test: Since the JM tests are somewhat elaborate we also used 
a simple null test for isotropy. Fig. [T] shows that the data naturally split 
into two separate groups: a group lying in the region 100 < RA < 250 (85 
sources, set 1), and another group (72 sources, set 2). For the test we fit the 
isotropic cosmo logical parameter Qm separately for each set, constraining 
n = + fiy = 1. The best fit for set 1 gives n^jS = 0.29 ± 0.05, while 
set 2 gives = 0.33 ± 0.05. This test yields no evidence for a systematic 
difference between the two sets. 

However it is reasonable to make a search over different directions defin- 
ing hemispheres labeled "up" and "down", fitting ^M-up and flM-down, in 
obvious notation. Such a search introduces two additional parameters, the 
orientation of defining "up" . This test finds a maximum difference 

^M-up — ^M-down = 0.14, with the axis position RA = 307°, Dec = —20°, 
roughly orthogonal to both celestial and galactic poles. Given that the one 
standard deviation error AQm = 0.05, and taking into account the 2 axis 
parameters, this difference is statistically significant at the 5 % level. That is, 
there is only 5 % probability that such a difference may arise by a fluctuation 
in an isotropic sample. 

While blunted by the Qm fltting procedures, the two- hemisphere test in- 
dicates a weak signal of anisotropy. This is quite consistent with its crudeness 
and lower sensitivity compared to the JM test. 

^Using the rescaled extinctions of Ref. J28] will not change the result 
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data set 




% 


Jo 


P 


X 


p 

r 


gold and silver - 


best fit 


-0.77 


1.96 


0.02 


228.4 


0.008 


gold and silver 




-0.70 


1.60 


= 


230.1 


0.010 


gold - best fit 




-0.68 


1.50 


0.021 


174.5 


0.092 


gold 




-0.63 


1.23 


= 


176.1 


0.11 


gold and silver - 


(5-best fit 


-0.53 


0.72 


0.029 


205.5 


0.085 


gold and silver - 


5 


-0.51 


0.57 


= 


209.2 


0.082 


gold - 6 -best fit 




-0.43 


0.28 


0.031 


152.4 


0.43 


gold - 6 




-0.45 


0.29 


= 


155.3 


0.43 



Table 2: The values for the gold and gold-silver set. Results are given for 
both the fit to the FRW model (p = 0) and the anisotropic Godel-Obukhov 
model. We also give results using the corrected extinction values, with the 
rescaling parameter 5 = —0.43. 

Model-Independent Correlation Summary To summarize results up 
to here: The supernova data are hardly consistent with isotropy when exam- 
ined with the model-independent angular correlation tests. At the same time 
there is at most a modest signal of fundamental anisotropy, which might be 
brought out further if the Godel-Obukhov model is appropriate. 

3.2 Fits 

Here we fit data to the Godel-Obukhov anisotropic cosmological model, Eqs. 
3.2.1 Nominal Extinctions 

In this study we use the extinction values published by Riess et al0. Pa- 
rameters are fit by minimizing x^, defined by 

^ t my ' 

where /Iq is the observed distance moduli, fip is the theoretical prediction 
(Eq. [6]) and (J/Xg are the reported errors. We vary six parameters: three FRW 
parameters fipo, deceleration (go) and jerk (jo), the anisotropy parameter 
p, and the two axis parameters. Results are summarized in Table [2l The 
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minimum in the rotating model for the gold and silver set is found to be 
228.4 for p = 0.02. Given six parameters this gives a x^/dof = 1.27, where 
dof stands for the degrees of freedom. This is to be compared with the 

= 230.1 in the FRW model, which has three parameters, and x^/dof = 
1.26. The rotating model provides a slightly worse overall fit to the data. 
Yet both fits are poor: the probability is quite low for a good fit to give 
such x^/dof values for so many data points. P- values listed in Tableware 
the probability to see the statistic or larger in x^, where v is the number of 
degrees of freedom. 

Dependence of x^ with p is shown in Fig. [2l From the fit we obtain a 
limit on the parameter p: one standard deviation limits are 

p < 0.046 (gold and silver); 

p < 0.049 (gold) . (10) 
3.2.2 Corrected Extinctions 

Jain and Ralston [28] found that re-scaling the host extinction values gener- 
ates a much better fit to the conventional isotropic model, decreasing the x^ 
values by more than 22 units, and bringing x^ /dof close to unity. Re-scaling 
also eliminates an alarming correlation between residuals and extinctions. 
The re-scaling correction replaces extinction values Ay by the rule 

Ay^Ay{5) = {l + 5)Ay, (11) 

with 5 = —0.43. Using the corrected Ay (6) values, the x^ values for the gold 
set (Table [2]) have x^ /dof = 1.01. These values are slightly different from 
that obtained in Ref. p^, since here we include the deceleration and jerk 
jo parameters, instead of Qm and and the constraint flu + = 1 is 
not imposed. Values for the gold and silver set with various assumptions are 
given in Table [2l 

With corrected extinctions the fit to angular dependence is just as good 
as the fit to the FRW metric. For the Gold data set the axis position is 
RA = 229.5° and Dec = —53.2°, roughly orthogonal to the galactic pole. The 
axis position for the Gold and Silver set is RA = 199.6° and Dec = —47.5°. 
We determine the upper limit of the parameter p by demanding that ^oi 
the anisotropic model is at most one unit larger than that in the FRW model, 
yielding 

p < 0.066 (gold and silver); 
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p < 0.067 (gold) 



(12) 



3.2.3 Quadrupole Model 

Up to here we have considered bounds on a specific model which incorporates 
rotation along with expansion. Yet the bound is really applicable to any 
model which at leading order gives a dipole contribution, A • f, where A is a 
preferred axis. 

A natural generalization of the model would incorporate a quadrupole 
contribution. We restrict our study to dependence going like cos^ a , while 
allowing the parameters of the preferred axis to vary. The anisotropic cor- 
rection to the distance modulus is then given by 

p' cos^a 

where p' can be positive or negative. In order to compare this with the 
distance modulus in the Godel-Obukhov model we may expand the last term 
in Eq. [6] to leading order in p. For the gold and silver set the minimum = 
225.1, x^/dof = 1.25 with parameters go = —0.68, jo = 1-5 and p' = 0.078. 
The fit obtained here is slightly better than that obtained in the FRW model. 
The previous bound becomes considerably looser, yielding —0.10 < p' < 0.17. 
If we use the extinction correction we find —0.14 < p' < 0.11. The minimum 
= 206.1, x^/dof = 1.145 with parameters go = —0.585, jo = 0.949 and 
p' = -0.060. 



3.3 Possible Selection Effects 

As seen in fig. 1 the data does not have a uniform angular distribution. Nor 
does one expect such a distribution given the selection effects of observations. 
The redshift values also show a correlation with angular position due to 
selection effects. We next ask the question how these effects might change 
our limits. 

We find that the selection effects do not affect our limits. Although 
selection effects might certainly generate an artificial signal of anisotropy, at 
least for some observables, the converse of hiding a true signal and showing 
a false null is a much different question. Much depends on how the analysis 
is arranged. In the present case these can generate neither a spurious signal 
nor a spurious null. 
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We demonstrate this with a simulation. We generated a simulated data 
set containing 400 sources. We arbitrarily assign an error of 0.2 to each sim- 
ulated distance modulus. The data is generated from a gaussian distribution 
with the mean value at any redshift given by the distance modulus obtained 
from the FRW metric. For 200 sources we uniformly assign a redshift range 
between and 1. For the remaining 200 sources we assign redshift values 
uniformly in the range to 2. These two sets are assigned two different, ar- 
bitrarily chosen, angular positions. This of course is an egrregiously selected 
data set. 

We fit the simulated data with our model for distance modulus seeking 
angular dependence. The analysis shows no significant signal of anisotropy. 
Despite the fact we have taken a highly anisotropic sample distribution on 
the sky, and a set that shows a very large angular redshift dependence, the 
correlation of distance modulus and angular position is unaffected. 

We carried out several examples of such simulations with different choice 
of angular positions and with data scattered at more than two positions on 
the sky. In no case do we find any significant signal of anisotropy. 

This shows that our tests of the correlation of the underlying distributions 
do not respond to substantial bias in the marginal distributions. In order for 
selection on the marginal distributions to conceal an actual correlation, one 
would have to deliberately cut out (select away) those special data capable 
of showing correlation. The only way to do so without cheating would be to 
select a very tiny range of redshifts so that the lever arm to find correlation 
would be too small. Such a cut will only make our limits more conservative. 
However they will still be rigourously applicable. 

4 Conclusions 

We have examined the magnitudes of large redshift type la supernova data 
for consistency with the Godel-Obukhov model. Using the nominal extinction 
coefficients gives poor fits, in which anisotropy makes the fit worse. Using 
corrected extinction coefficients gives excellent fits, while also opening up the 
allowed region of the anisotropy parameter. We have also tested for dipole 
and quadrupole anisotropy in the supernova data. In none of the cases have 
we found strong evidence in favor of any anisotropic effects. It is nevertheless 
impossible to rule out anisotropy entirely. 
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Figure 2: for the 186 object data set for the Godel-Obukhov metric as a 
function of the parameter p. 
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